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ABSTRACT 



Context. Infrared dark clouds (IRDCs) likely represent very early stages of high-mass star/star cluster formation. 
Aims. The aim is to determine the physical properties and spatial distribution of dense clumps in the IRDC MSXDC G304.74+01.32 
(G304.74), and bring these characteristics into relation to theories concerning the origin of IRDCs and their fragmentation into clumps 
and star-forming cores. 

Methods. G304.74 was mapped in the 870 yum dust continuum with the LABOCA bolometer on APEX. The 870 yum map was 
compared with the 1.2 mm continuum map of the cloud by Beltran et al. (2006). Archival MSX and IRAS infrared data were used to 
study the nature and properties of the submillimetre clumps within the cloud. The H2 column densities were estimated using both the 
870 yum dust emission and the MSX 8 yum extinction data. The obtained values were compared with near-infrared extinction which 
could be estimated along a few lines of sight. We compared the clump masses and their spatial distribution in G304.74 with those in 
several other recently studied IRDCs. 

Results. Twelve clumps were identified from the 870 yum dust continuum map. Three of them are associated with the MSX and 
IRAS point sources. Moreover, one of the clumps (SMM 6) is associated with two MSX 8 yum point-like sources. Thus, there are 8 
clumps within G304.74 which are not associated with mid-infrared (MIR) emission. The H2 column densities derived from the dust 
continuum and extinction data are similar. The comparison suggests that the dust temperature may be elevated (20-30 K) near the 
southern end of the cloud, whereas the starless clumps in the centre and in the north are cool (T^ ~ 15 K). There is a high likelihood 
that the clump mass distributions in G304.74 and in several other IRDCs represent the samples of the same parent distribution. In 
most cases the spatial distributions of clumps in IRDCs do not deviate significantly from random distributions. 
Conclusions. G304.74 contains several massive clumps that are not associated with MIR emission. On statistical grounds it is likely 
that some of them are or harbour high-mass starless cores (HMSCs). The fact that the clump mass distributions (resembling the high- 
mass stellar IMF), and in some cases also the random-like spatial distributions, seem to be comparable between different IRDCs, is 
consistent with the idea that the origin of IRDCs, and their further sub-fragmentation down to scales of clumps is caused by supersonic 
turbulence in accordance with results from giant molecular clouds. 

Key words. Stars: formation - ISM: clouds - ISM: structure - Radio continuum: ISM - Submillimetre 



1. Introduction 

1.1. Infrared dark clouds 

The so-called infrared dark clouds (IRDCs) were discovered by 
the Infrared Space Observatory (ISO; IPerault et al. 1996] l and 
the Midcourse Space Experiment (MSX; Eganetal. 1998 1; the 
clouds were detected as dark absorption objects agamst the 
bright mid-infrared (MIR) Galactic background radiation. Based 
on the MSX 8 fim data of the Galactic plane from I - - 360° 
and \b\ < 5°, Simon et al. (2006a) identified almost 1 1 000 IRDC 
candidates. Recently, Peretto & Fuller (2009) used Spitzer satel- 
lite data to extract about 9000 new IRDC candidates. 

Molecular line and dust continuum studies of IRDCs have 
shown that they are cold (T < 25 K), dense (n(H2) > 10^ 
cm-^ A^(H2) > 10^2 cm-2), and massive ( ~ 10^ - 10^ Mq) 
structures with sizes of ~ 1 - 15 pc (e.g., |Carey et al. 1998[ 
ISimon et al. 2006bl IRathborne et al. 20061 hereafter RJS06; 



Du & Yang 20081 |Vasyunina et al. 2009|. Most IRDCs are fila- 
mentary (e.g., IPeretto & Fuller 20 09), and contain density en- 
hancements, or clumpfl that are visible in (sub)millimetre 
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dust continuum m aps (e.g., Carey et al. 2000 Garay et al. 2004 
lOrmel et al. 20051 RJS06). The cold clumps (i.e., clumps unas- 
sociated with MSX 8 /im emission) identified by RJS06 have 
typical sizes and masses of ~ 0.5 pc and ~ 120 M0, respec- 
tively. Because IRDCs have clumpy structures, they are likely 
to be in an early stage of fragmentation (e.g., RBG09 and refer- 
ences therein). 

The radial galactocentric distribution of IRDCs peaks at 
Rcc = 5 kpc in the 1st Galactic quadrant, and at Rqc - 6 
kpc in the 4th quadrant, which correspond to the location of 
the Scutum-Centaurus spiral arm (see ISimon et al. 200661 and 
IJackson et ar2 008). This, together with the fact that IRDCs 
have sizes and masses similar to those of warm massive cluster- 
forming regions (e.g., ILada & Lada 20031 IMotte et al. 2003l l. 
has led to the suggestion that IRDCs represent the very 



' We prefer to use the term "clump" according to e.g., Ragan et al. 
(2009, hereafter RBG09), rather than the term "core" which was used 
by e.g., RJS06. Clumps have masses and sizes (radii) of ~ 10- 10^ Mq 
and ~ 0.1 - 1 pc, respectively (e.g., |Bergin & Tafalla 2007| l. 
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early stages of high-mass star/star cluster formation. Indeed, 
several studies have found signs of ongoing star formation 
within IRDCs. These include CH3OH and H2O masers 
( Beuther et al . 2002bl IPillai et al. 2006al |Wang et al. 2006 
Ellingse n 2006|, outflow signatures dBeuther et al. 200? 
'Beuther & Sridharan 20071 ISakai et al. 2008b, and as- 
sociated infrared sources (e.g., IRathborne et al. 20051 
[Beuthe r & Steinacker 20071 [Chambers et al. 2009; RBG09). 
Cold, dense clumps are suggested to host/represent precursors 
of hot molecular cores, i.e., high-mass starless cores (HMSCs; 
e.g., ISridharan et al. 20051 IBeuther et al. 20071). A few hot 
cores have been found within IRDCs (Rath borne et al. 20071 
2008), but it is important to increase the sample in order 
to study the sequence of events. It should be noted that the 
formation of hot cores may not be universal in IRDCs - some of 
them seem to form only low- to intermediate-mass stars (e.g., 
[van der Wiel & Shipman 2008| . 

1.2. Infrared dark cloud MSXDC G304.74+1.32 

The IRDC studied in this paper was designated as MSXDC 
G304.74+1.32 (hereafter, G304.74) by Simon et al. (2006a). 
G304.74 was observed by Beltran et al. (2006) with the SIMBA 
bolometer array on SEST at 1.2 mm. Beltran et al. (2006) identi- 
fied eight millimetre clumps within the cloud, out of which four 
were found to be associated with MSX point sources. The dis- 
tance to the cloud, 2.4 kpc, is a kinematic distance estimated 
from the velocity of the CS line observed by Fontani et al. 
(2005). The galactocentric distance is ~ 7.4 kpc. 

G304.74 was chosen for the present study because its rela- 
tively close distance allows for reasonably good spatial resolu- 
tion in order to study its substructure. Moreover, a relatively high 
number of clumps (8) were already identified from the cloud 
(see above), and thus it was considered an appropriate object for 
studying IRDC fragmentation. 



tors for LABOCA), and the star CW Leo as secondary calibra- 
tor. The uncertainty due to flux calibration was estimated to be 
~ 10%. The telescope focus and pointing were checked using 
the planet Saturn, the star rj Carinae, and the H2O maser source 
B13134 (305.80-00.24). The submm zenith opacity was deter- 
mined using the sky-dip method and the values varied from 0.20 
to 0.38, with a median value of 0.30. The observations were thus 
conducted in fair weather conditions. 

The observations were made using the the on-the-fly (OTP) 
mapping mode with a scanning speed of 3' s The size of the 
OTP map is 34' x 34'. This was achieved by 204 subscans of 34' 
in length (parallel with the R.A. axis), spaced by 10". The area 
was mapped five times, with a total observing time of 5.2 h. In 
this manner, a uniform rms noise level of cr = 0.03 Jy beam ' 
was reached in the central ~ 15' x 15' area which covers the 
target source. 

The data reduction was performed using the BoA (Bolometer 
Array Analysis Software) software package according to guide- 
lines in the BoA User and Reference Manual (2007j5 The data 
reduction included flat-fielding, flagging bad/dark channels and 
data according to telescope speed and acceleration, correcting 
for the atmospheric opacity, division into subscans, baseline sub- 
tractions and median-noise removal (correction for sky noise), 
despiking, and filtering-out of the low frequencies of the l/f- 
noise. Pinally, the five individual maps were coadded. 

We note that even the source model was iteratively used in 
the reduction process, the resulting final map (Pig.[T] left panel) 
has negative artefacts ("holes") around regions of bright emis- 
sion. The depths of these negative holes are ~ 15 - 50% (~ 30% 
on average) of the nearest peak object brightness, and thus they 
are likely to introduce additional uncertainty in the source flux 
densities. We estimate that the total flux density uncertainty due 
to calibration and negative artefacts is in the range ~ 20 - 50%. 
These uncertainties are not, however, taken into account in the 
analyses presented in this paper 



In this paper, we present the results of our 870 yum dust con- 
tinuum mapping of G304.74. The paper pays special attention to 
the clumpy structure of the cloud, and the clump mass distribu- 
tion, and thereby addresses the fragmentation of IRDCs. The ob- 
servations and data-reduction procedures are described in Sect. 
2. The observational results are presented in Sect. 3. The MSX 
8 jum optical thicknesses toward the submm peaks are derived in 
Sect. 4. In Sect. 5, we describe the methods used to derive the 
physical properties of the observed clumps. In Sect. 6, we dis- 
cuss the results of our study, and in Sect. 7, we summarise our 
main conclusions. 



2. Observations and data reduction 

2.1. Submillimetre dust continuum 

The 870 yum dust continuum observations toward G304.74 were 
carried out on 25-27 April 2009 with the 295 channel bolometer 
array LABOCA (Large APEX Bolometer Camera) on APEX. 
The central frequency of the instrument is 345 GHz, and the 
bandwidth is 60 GHz. The half -power beam width (HPBW) of 
the telescope is 18'.'6 (~ 0.22 pc at 2.4 kpc) at the frequency 
used. The total field of view (PoV) of LABOCA is 1 1 '.4. The in- 
strument and its observing modes are described in Siringo et al. 
(2009). 

Absolute flux calibration was achieved through observations 
of the planets Mars, Uranus, and Neptune (the primary calibra- 



2.2. Complementary near-, mid-, and far-infrared data 

The target source (I = 304°74, b = 1°32) is not included in the 
GLIMPSE (3.6, 4.5, 5.8, 8.0 /zm) and MIPSGAL (24 and 70 fim) 
surveys of Spitzer which cover the Galactic latitudes \b\ < 1°. 
We have used the near-infrared 2MASS (J, H, and K^) sur- 
vey data archive (Skrutskie et al. 2006)1, data products from the 
MSX survey (Price et al. 2001: Egan et al. 2003| ), and the IRAS 
(Infrared Astronomical Satellite) MIR and far-infrared (PIR) ob- 
servations. 

3. Observational results 

3. 1. Clump identification 

The obtained LABOCA map is presented in the left panel of 
Pig. [1] The cloud has a filamentary appearance and extends 
over about 13'.2, i.e., about 9.2 pc at the cloud's distance. In or- 
der to identify clumps in the LABOCA map, we employed the 
two-dimensional clumpfind algorithm, cl£ind2d, developed 
by Williams et al. (1994). The cl£ind2d routine determines the 
peak position, the PWHM (full width at half maximum) size (not 
corrected for beam size), and the peak and total integrated flux 
density of clump based on specified contour levels. The algo- 
rithm requires two input parameters: 1) the intensity threshold, 
i.e., the lowest contour level, which determines the minimum 

^ |http : //www ■ astro . uni-bonn . de/boawiki/Boa] 
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emission to be included into clumps; and 2) the stepsize which 
determines the required "contrast" between two clumps to be 
considered as different objects. We set both the intensity thresh- 
old and the stepsize to 2cr of the background noise level. With 
these settings, G304.74 divides into 12 clumps. 

We note that the total number and size of clumps identified 
by clumpfind is quite sensitive to the selected contour levels. 
The selected 2cr contour levels turned out to give the best agree- 
ment with the identification by eye. This selection is also recom- 
mended by WilHams et al. (1994). See also Pineda et al. (2009) 
for a recent discussion of clump identification with clumpfind. 
Pineda et al. (2009) concluded that small changes in the thresh- 
old and/or stepsize values can lead to important changes in the 
number of identified clumps. 

3.2. Observed properties of the clumps 

All the clumps have a peak flux density > 5cr (i.e., > 0.15 Jy 
beam ') relative to the local background. The coordinates, peak 
and integrated flux densities, and deconvolved angular FWHM 
diameters (^s) are listed in Cols. (2)-(6) of Table [T] The values 
of 6s correspond to the beam corrected FWHM size. Column (7) 
lists the eff'ective radius, 7?efF- This radius is defined as /?eif = 
y/A/n, where A is the projected area with the clump boundaries. 
The corresponding circles are overlaid on the LABOCA map 
in the right panel of Fig. [1] They are centred on the dust peak 
positions (Cols. (2) and (3)). Note that the clumpfind algorithm 
makes no assumption about the shape of the clump. 

The left panel of Fig. |2] shows the wide-field MSX 8 yum 
image (extracted from the MSX Galactic Plane Survey imaged 
towards G304.74. The MSX 8 fim image zooming to the dark 
cloud with LABOCA contours is shown in the right panel of 
Fig.|2] There is a good agreement between the morphologies of 
the 870 yum continuum emission and the 8 /im extinction visible 
as a dark lane in the MSX image. The resolution of the MSX 
image, 18'.'3, is similar to that of the 870 fim LABOCA map. 

Three of the submm clumps are associated with sources both 
from the MSX and IRAS point source catalogues (within ~ 6 - 
23" and ~ 4- 20" from the dust peak position, respectively; see 
Fig.[T] left). These are designated according to the IRAS names. 
The remaining nine clumps are named SMM 1, SMM 2, etc. One 
of the submm clumps, SMM 6, appears to be associated with two 
8 fim sources from the MSX catalogue. The stronger 8 fim source 
is located at ~ 6" from the dust peak position (0.46 Jy beam"'), 
whereas the weaker one is located at ~ 19" from the closest dust 
emission peak (0.41 Jy beam"') within SMM 6. Furthermore, 
the eastern edge of SMM 4 is bright in 8 fim as the MSX source 
associated with IRAS 13037-61 12 extends about 55"west. Thus 
there are seven clumps which are completely dark in the MIR. 

The flux densities at MIR and FIR wavelengths retrieved 
from the MSX (8.28, 12.13, 14.65, 21.34 ^/m) and IRAS (12, 25, 
60, 100 fim) archives for all the IRAS sources and for the two 8 
fim point sources associated with SMM 6 are listed in Tables |2] 
and |3] respectively. 



4. IRDC extinction and the 8 fim optical thicknesses 

4.1. Principle 

can be 



Table 3. IRAS point sources in the IRDC G304.74. 



The observed 8 fim intensity towards an IRDC, I\ 
written as (e.g.. ,Bacmann et al. 2000) 



obs 

IRDC 





5 12 


525 


S 60 


S 100 


Name 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


IRAS 13037-6112 


1.84 


7.41 


64.96 


< 196.5 


IRAS 13039-6108 


3.99 


7.43 


105.6 


196.5 


IRAS 13042-6105 


< 1.31 


0.67 


< 4.97 


< 196.5 



robs 

'IRDC 



Ibi,e-'""" + If, 



(1) 



where /bg is the 8 fim intensity of the background emission, Tg^m 
is the 8 fim optical thickness, and /foie is the 8 fim intensity con- 
tribution from foreground material. Equation ([TJ can be inverted 
to get the value of Tg^n, as 



Tg/jm 



•In 



'IRDC 



(2) 



Following the notation of Bacmann et al. (2000) and Peretto 
& Fuller (2009), the observed MIR intensity around the IRDC is 



I MIR 



ho + Ifoi 



(3) 



To obtain an estimate for the value of /mir, we used the technique 
of spatial median filtering jSimon et al. 2006al ). In this method, 
^MiR for each pixel is estimated by taking the median of intensi- 
ties inside a surrounding filter region. The size of the filter must 
be larger than the size of IRDC; the method will only capture 
emission fluctuations on scales larger than the size of the cloud. 
Given the size of G304.74 (> 13', i.e., over 130 6"pixels of the 
MSX map), we used a circular spatial filter with a radius of 15' 
(the same as used by Simon et al. (2006a)). We note that the to- 
tal size of the image used in the filtering process was 0°5 x 0°5 
(see Fig. |2] left). Figure [3] shows the obtained model of /mir. 
Both the mean and median value of /mir are 1.23 + 0.02 x 10"^ 
W m"^ sr"', where the uncertainty is the standard deviation of 
the pixel values. The above value was adopted in the following 
analysis. 



ZOmOO.Os 



22m 00.0$ 



24m 00.0s 



26m 00.0s 



-61d30m00.0$ 



30.0s 13h07m00.0s 06m30.0s 
Right ascension (J2000) 

Fig. 3. Median filter estimate of /mir - The intensity scale is linear 
from 1.20 X lO"*" to 1.26 x lO"*" W m"^ sr'', wifli both mean and 
median being 1.23 + 0.02 x 10"*' W m^^ sr"'. 
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Fig. 1. Left: LABOCA map of the 870 fim dust continuum emission from the IRDC G304.74. The green plus signs and white boxes 
mark the positions of the MSX 8 fim sources in the field (cf. Fig. |2]i and IRAS point sources, respectively. The intensity range is 
from -0.17 to 0.66 Jy beam as indicated in the colour bar. The beam HPBW (18'.'6) and scale-bar are shown in the bottom left. 
Right: Grey-scale image of 870 fj.m emission shown in the left panel with the locations and sizes of the submm clumps determined 
by the clumpf ind algorithm. The circles represent the area associated with each clump but not the shape. 




08m30.0s 13h07m00.0s 05ni30.0s sons i3hD7mDD.Ds D6m3D.Ds 



Right ascension (J2000) Right ascension (J2000) 

Fig. 2. Left: Wide-field (0^5 xO?5, i.e., ~ 21 x21 pc at 2.4 kpc) MSX 8 /jm image showing the MIR extinction of the IRDC G304.74. 
Riglit: Blow-up of the left panel showing the same region as in Fig.[T] The MSX 8 fim image is overlaid with contours of LABOCA 
870 fim emission. The intensity range of a grey-scale image is from 1.86 x 10"^ to 2.32 x 10"'' W m"- sr"'. The contours are 0.06 
(2cr) to 0.66 Jy beam in steps of 0.06 Jy beam '. The beam HPBW (18'.'3) and scale-bar are shown in the bottom left. 



4.2. Background and foreground estimation 

To estimate the contribution of /bg and /fore to the value of 
^MiR (Eq. (O, and Fig. [3J, we used the similar analysis as used 
by Butler & Tan (2009, their section 3.1). We assumed that 
the radial distribution of Galactic 8 fim PAH (polycyclic aro- 
matic hydrocarbon) emission has the same exponential form 
as the radial distribution of molecular gas surface density in 
the Galaxy (the radial scale length is Hr - 3.5 ± 1.0 kpc; 
IWilhams & McKee 19971 1 . We note that this distribution is sim- 
ilar to that of the Galactic surface density of OB associations 
(IMcKee & Willi ams 1997l l. which was used by Butler & Tan 
(2009, their Eq. (3)). The assumption that the radial distribu- 
tion of PAHs follows the molecular gas is supported by the re- 



sults of Mattila et al. (1999) who found that the UIR (unidenti- 
fied infrared) emission band intensities in the disk of the galaxy 
NGC 891 closely follow those of the CO emission and 1.3 mm 
dust continuum emission. However, UIR band intensities were 
found to be quite different from that of the Ha emission (see 
Fig. 5 of Mattila et al. (1999)). Moreover, Kahanpaa et al. (2003) 
found strong correlation between the Galactic UIR bands and 
CO emission. To estimate the vertical distribution of PAH emis- 
sion, we used the Galactic vertical distribution of CO emission 
and adopted the scale length z\/2 = 70 + 10 pc JDame et al. 1987t 
IBronfman et al. 1988HMaihotra 1994l l. 

We first calculated the so-called "foreground intensity ratio", 
/fore, which is defined as the ratio of the column of PAH emission 
between the observer and the cloud and the total column through 
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Table 1. Submillimetre clumps identified by the clumpfind algorithm in the IRDC G304.74. 



5 



Name 


Peak position 
0-2000.0 [h:m:s] ^2000.0 [°-'-"] 


870 — 

[Jy beam"'] 


870 

[ly] 


e, 
["] 


^cftl 1.2 mm clump No/_ 
["] (Beltran et al. 2006") 


SMM 1 


1 3 06 20 5 


-61 30 1 5 


0.46 


2.0 


36 


36 




SMM2 


1 3 06 26 9 


-61 29 39 


0.30 


1.0 


38 


28 




SMM 3 


13 06 35.8 


-61 28 53 


0.50 


3.4 


45 


42 


1 


SMM 4 


13 06 44.7 


-61 28 35 


0.66 


4.1 


38 


38 


3 


IRAS 13037-6112 


13 06 49.8 


-61 28 07 


0.47 


1.6 


24 


26 


2 


SMM 5 


13 06 49.9 


-61 26 55 


0.19 


0.7 


28 


26 




SMM 6 


13 06 51.1 


-61 27 49 


0.46 


1.8 


27 


30 


4 


SMM 7 


13 07 02.6 


-61 26 18 


0.26 


0.9 


24 


30 


8 


IRAS 13039-6108 


13 07 06.4 


-61 24 29 


0.47 


2.6 


37 


40 


5 


SMM 8 


13 07 08.9 


-61 23 25 


0.31 


1.1 


26 


29 




SMM 9 


13 07 12.7 


-61 22 49 


0.38 


1.6 


35 


32 


6 


IRAS 13042-6105 


13 07 20.3 


-61 21 45 


0.21 


0.8 


28 


27 


7 



" The peak and total flux densities have estimated uncertainties of ~ 20 - 50%, depending on the intensity of negative holes around the clump 
(see Sect. 2.1). 

* The effective radius, R^ff = s/A/n, where A is the area assigned to the clump, i.e., area contained within a 2tr contour The radius has not 
been deconvolved from the beam size. 

Corresponding 1.2 mm clump in Beltran et al. (2006). 



Table 2. MSX point sources in the IRDC G304.74. 







S 8.28 


S 12.13 


S 14.65 


S 21.34 


Associated clump 


MSX6C designation 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


IRAS 13037-6112 


G304.7728+01.3431 


0.876 ± 0.037 


< 0.876 ± 0.076 


1.027 ±0.076 


< 1.969 ±0.142 


SMM 6 


G304.7738+01.3522 


0.213 ±0.011 




< 0.753 ±0.061 


2.972 ±0.193 


SMM 6 


G304.7800+01.3597 


0.184 ±0.010 








IRAS 13039-6108 


G304.8074+01.4037 


1.875 ±0.077 


2.125 ±0.123 


< 0.730 ± 0.063 


< 1.678 ±0.129 


IRAS 13042-6105 


G304.8366+01.4472 


0.147 ±0.010 









the Galaxy along the same line of sight, up to the galactocentric 
distance 16 kpc. For G304.74, we obtained /fore = 0.514 ± 0.05, 
where the uncertainty represents the mimimum-maximum error 
based on the uncertainties in the scale lengths (see above). The 
distance of the Sun from the Galactic Centre was assumed to 
be 8.5 kpc in this calculation. The intensity of the background 
emission is then given by 

4g = (1 - /foi-e)^MIR • (4) 

The resulting values are hg = 5.98 + 0.62 x 10"^ W nr^ sr ' and 
^foie = 6.32 ± 0.65 X 10"^ W m"^ sr"'. The errors in /bg and /fore 
are propagated from the standard deviation of /mir and the un- 
certainty in /fore. The above value of /fore is in good agreement 
with the lowest observed 8 /urn intensities as demonstrated in 
Fig.|4] We note that the anti-correlation between the submm and 
MIR intensities allows for a determination of /bg and /fore (e.g., 
[Johnstone et al. 2003 Fig. 3 therein). However, the correlation 
shown in Fig.|4]is not very clear and thus only the /fore-method 
is considered in the present paper (see Sect. 6.1 for further dis- 
cussion). 

The observed 8 /im intensities toward the submm dust con- 
tinuum peak positions, and the 8 /jm optical thicknesses calcu- 
lated from Eq. (|2]l, are listed in Cols. (2) and (3) of Table|4] The 
formal error in Tg|^^^ was calculated by propagating the errors in 

/fore and /bg. In order to determine I°^qq values in the submm 
peak positions, the MSX image was gridded with 9'.'3 pixels, 
i.e., the LABOCA pixel size. We also estimated the possibil- 
ity that part of the observed 8 //m intensity originates from the 
bright surroundings of a MIR dark clump. The measured point 
spread function (PSF) in the MSX image is a Gaussian with a 



FWHM of 22". Convolution with this PSF reduces the breadths 
of dark filaments, and, for the narrowest of them, increases the 
minimum intensity in the middle. We estimate that for a com- 
pletely opaque source with a diameter of < 30", the MIR radia- 
tion "leakage" from the surroundings amounts to > 20% of the 
background level. This contribution was, however, neglected in 
the above analysis which uses only the peak (minimum) 8 jum 
intensities. Also, most of the MIR dark clumps are larger than 
30" (see Tableffl Cols. (6) and (7)). 

4.3. Dust temperature estimates 

Because dust continuum emission is optically thin at (sub)mm 
wavelengths, the 870 fim radiation intensity is given by 

/gvo - B870(rd)T870 , (5) 

where Bfi,jo{Td) is the Planck function with dust temperature T^. 
According to the Ossenkopf & Henning (1994, hereafter OH94) 
dust model used in the present paper (see Sect. 5.1 for more de- 
tails), the ratio of dust opacities per unit dust mass at 8.^ and 
870 fim, /fg/zcgvo, and thus the correponding ratio of the optical 
thicknesses, r8/T87o, is about 865. The obtained ratio is consis- 
tent with the results of Johnstone et al. (2003) and Ormel et al. 
(2005), who found that (assuming = 15 K) /cg/fsso ~ 640 
and 870, respectively. By using this ratio, and the rg values de- 
rived above, estimates for the dust temperatures, Td, towards the 
submm peaks can be derived from Eq.(|5]l. The resulting val- 



Note, for the MSX band A filter the 50% cutofi's are at 6.8 and 
10.8 jim (central wavelength is 8.8 fim), whereas the isophotal central 
wavelength of the filter is 8.28 fim. 
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Table 4. The observed 8 jum intensity, the peak optical thickness, 
and the dust temperature toward the submm peak positions. 



Name 


robs 
^IRDC 

[10-' sr'] 




Ti 
[K] 


SMM 1 


9.89 


0.52 ±0.21 


30.6 ± 17.0 


SMM2 


10.39 


0.38 ±0.19 


28.0 ± 19.6 


SMM 3 


11.73 


0.10^ 




SMM 4 


10.97 


0.25 ±0.17 


J 


SMM 5 


8.03 


1.25 ±0.40 


10.1 ±2.8 


SMM 7 


8.08 


1.22 ±0.39 


12.0 ± 3.7 


SMM 8 


7.58 


1.56 ±0.53 


11.5 ± 3.7 


SMM 9 


6.99 


2.19 ±0.98 


10.8 ±4.3 



" No error given because it is greater than the actual value. 
* The value of could not be reasonably estimated. 



3x1 0" 



2x10" 



0.0 



0.2 



0.4 
(Jy beam" 



].6 



0.8 



ues are Hsted in Col. (4) or Table |4] The quoted errors are the 
minimum-maximum errors derived from the uncertainty in Tg. 

The uncertainties of the estimates are very large in the 
southern part (SMM 1, SMM 2; in the case of SMM 3 and 
SMM 4 no reasonable estimate could be obtained because of 
the large error in rg). On the other hand, the estimates with mod- 
erate formal errors in the central and northern part of the cloud 
suggest low temperatures of slightly above 10 K. It should be 
noted, however, that there are large uncertainties concerning the 
dust opacities and the determination of the rg. This issue will be 
raised again in Sect. 6.2. 

In view of these uncertainties, dust temperatures estimated 
above should be compared with previous temperature deter- 
minations from molecular lines. Observations towards sev- 



eral other IRDC cl umps (H2CO, [Carey et al. 1998| CH3CCH 
Teyssier et al. 2002 ; NH^. Pillai et al. 2006bl ISakai et al. 20"08] l 
provide gas kinetic temperatures of Jkin ~ 10 - 20 K. The 
assumption that - Tkin is probably valid in dense clouds 



(n(H2) > 10^ cm"^; e.g.. Goldsmith & Langer 1978). Based on 



these considerations, we assume in what follows that in gen- 
eral Td = 15 K. Exceptions are made in the cases of two IRAS 
sources, for which temperatures can be derived from the spectral 
energy distributions (SEDs). We will also separately consider the 
possibility that is elevated in the four southern clumps SMM 
1-4. 



5. Physical parameters of the clumps 

5.1. Spectral energy distributions of tine IRAS sources 

The MSX, IRAS, LABOCA, and SIMBA data were used to 
fit the spectral energy distributions (SEDs) of IRAS 13037- 
6112 and IRAS 13039-6108. The SIMBA 1.2 mm flux den- 
sities of these sources are 0.89 and 1.36 Jy, respectively 
dBeltran et al. 20061 Table 2 therein). Note that there is not 
enough data points for IRAS 13042-6105 in order to con- 
struct a reasonable SED (e.g., most of its IRAS flux den- 
sities are only upper limits, see Table O. The derived 
SEDs are shown in Fig. |5] The least-squares fitting rou- 
tine used in the derivation of the SEES minimises the sum 

i:f=i [logi()(5f^(^/)) - logio(5r'(/l/))]^ where is the num- 
ber of data points (7 and 8 for IRAS 13037-6112 and 13039- 
6108, respectively), Sf"^(Ai) is the observed flux density, and 
S ™'^(/l,) is the model flux density. In both cases, the data were 
fitted by a two-temperature composite model. It was assumed 



Fig. 4. MSX 8 iim versus LABOCA 870 //m intensity for 
G304.74. The horizontal line marks the 8 fim foreground inten- 
sity, /fore - 6.32 X 10"^ W nr^ sr"', derived from the foreground 
intensity ratio, /fore (see main text). 



that both components at different temperatures emit as a black- 
body modified by the wavelength-dependent dust opacity, (see 
below). The best-fit model SEDs overestimate the flux densities 
at ~ 12-20 jum, but underestimate them at ~ 8 and 25 fim. It 
should be noted that the flux densities included in the SEDs are 
measured using telescopes with different beam sizes. Thus the 
flux densities obtained for extended sources are not fully com- 
parable, and this can in part explain discrepancies between MSX 
(18'.'3) and IRAS (~ 2' at 12 fim to ~ 4' at 100 fim) fluxes at 12 
and - 20- 25 fim. On the other hand, SIMBA and LABOCA 
flux densities at A - 1.2 mm and 0.87 mm refer to clump ar- 
eas (Rgff is typically ~ 30") derived by cltompfind. /?eif values 
are similar for both IRAS 13037-6122 and IRAS 13039-6108 
(~ 30 - 40" or ~ 0.35 - 0.47 pc). Assuming that the emission in 
the IRAS bands is confined in the region of the submm clump, 
the characteristic spatial scale associated with the SEDs is < 0.5 
pc. 

We have adopted a dust-to-gas mass ratio of = 1/100, 
a value which has often been used in the IRDC studies (e.g., 
RJS06; ' Vasyunina et al. 2009t IParsons et al. 2009l l. However, 
this value can differ from 1/100. For instance, Draine et al. 
(2007) determined a value of ^ 1/186 based on observed de- 
pletions of heavy elements in the Galaxy. Dust opacities we have 
adopted correspond to a MRN size distribution with thick ice 
mantle^ at a gas density of «h = 10^ cm"-' (OH94). The result- 
ing SED parameters are given in Table |5] The total (cold+hot) 
mass and the integrated bolometric luminosity are given in 
Cols. (2) and (3) of Table|5] respectively. The temperatures of the 
two components are listed in Cols. (4) and (5). In Cols. (6) and 
(7), we give the mass and luminosity fractions of the cold com- 
ponent versus the total mass and luminosity. Column (8) lists 
the mass-to-luminosity ratio, Mtot/iboi, which is an evolutionary 
indicator of the clump as it is expected to decrease with time. 



The SED fitting routine was originally written by J. Steinacker. 



* In cold, dense interiors of the IRDCs, dust grains are supposed to 
be coagulated and covered by icy mantles (e.g., Peretto & Fuller 20091 
and references therein). This is supported by e.g., the observed 
depl etion of H2C O ( [Care y et al. 19 98, 2000), CO (JPillai et al. 2007] 
Zha ng et al. 2009^ , and CS (Beut her & Henn ing 20091 in IRDCs. 
Moreover, Butler & Tan (2009) found some evidence for dust opacity 
changes within IRDCs, which could be caused by ice mantle formation 
and grain growth. 
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The envelope mass decreases during the star formation process, 
and the luminosity of the embedded star or stellar cluster rises 
(e.g.,]Sndharan et al. 2002). We note that the adopted dust opac- 
ity model (corresponding to particles coated with thick ice man- 
tles) is not likely to be appropriate for hot dust, and therefore 
the total bolometric luminosity for the hot component should 
be taken with caution. The MIR spectral features (such as PAH 
emission) also cause the fit to the hot component being more un- 
certain than the fit to the cold part of the spectrum. However, this 
should not alter the fact that for both IRAS sources the mass of 
the hot component is negligible (~ 1 - 2 x 10"^ Mq), and thus 
the bulk of the material is cold (Mcoid/Mtot ~ 1). The bolomet- 
ric temperature of IRAS 13039-6108 (T^oW - 22 K) is in good 
agreement with the rotation temperature of 18 K derived from 
C'^O by Fontani et al. (2005). 



IRAS 13037-6112 



60 M 




10000 



IRAS 13039-6108 



10" 



™ 10"'"- 



T -?'^ K 

L„„„ = l.(ixl0' L 



10" 



10" 




100 
X [/im] 



1000 



10000 



Fig. 5. Best-fit SEDs of IRAS 13037-6112 and IRAS 13039- 
6108. MSX data points are shown by circles, IRAS data points 
are shown by diamonds, and LABOCA and SIMBA (sub)mm 
data points are indicated by squares. Upper limits are indicated 
by arrows. The solid lines in both plots represent the sum of two 
(coldH-hot) components (see Cols. (4) and (5) of Table |5]l. The 
absorption features appearing at ~ 3.1 and 10 are attributed 
to the H2O ice and silicates, respectively (OH94 and references 
therein). 



5.2. Linear sizes, masses, and H2 number densities 

The clump linear radii were computed from the effective radii 
listed in Col. (7) of Table[T] The clump masses (gas-i-dust mass. 
Meant) were calculated from their integrated 870 yum flux density, 
5870 (Table[T] Col. (5)), assuming that the thermal dust emission 
is optically thin ( |Hildebrand^l983) : 



S 870^/^ 



(6) 



where d is the source distance. For IRAS 13037-6112 and 
13039-6108, we assumed the dust temperatures to be the same 
as their bolometric temperatures, Tcoid - 22 K, resulting from 
SED fits (see Table|5] Col. (4)). For all the other clumps (SMM 
1, SMM 2,..., and IRAS 13042-6105), it was assumed that 
Td = 15 K (Sect. 4.3). We assumed that k^) - 0.17 m^ kg"'. 
This value is interpolated from OH94, see Sect. 5.1. As in the 
SED fits, the value 1/100 is adopted for R^. 

The volume-averaged H2 number densities, («(H2)), were 
calculated assuming a spherical geometry for the clumps, using 
the formula 



<n(H2)) : 



(P) 



(7) 



where (p) = MconJ {^/^^^Ift) '^^e matter density, yUe, - 2.8 
is the mean molecular weight per H2 molecule (assuming a 
10% helium abundance), and me is the mass of the hydrogen 
atom. The obtained radii, masses, and volume-averaged H2 num- 
ber densities are given in Cols. (2), (3), and (7) of Table |6] 
respectively. The typical density in the clumps is likely to be 
higher than the volume-averaged value because of substructure. 
Note that our masses are in general larger than those derived by 
Beltran et al. (2006) who assumed a dust temperature of 30 K 
and a /(■i.2mm of 0.1 m^ kg"' in their mass estimates. They also 
obtained higher densities because of using the FWHM radii (in- 
stead of effective radii) and a mean molecular weight of 2.29. 

5.3. H2 column densities 

The H2 column densities, A^(H2), towards the submm peaks were 
estimated using 1) the dust emission at 870 fim, and 2) the dust 
extinction at 8 yum. In the first method, the LABOCA intensities 
were converted to A^(H2) using the equation 



where 7^"^ = ^g^jJ'/Qbeam is the observed dust peak surface 
brightness (fibeam is the beam solid angle), which is related to 
the peak flux density via 1 Jy/18'.'6 beam = 1 .085 x 10"'^ W m"^ 
Hz"' sr"'. The dust temperature values used were identical to 
those adopted in the mass estimates (Eq. (|6]l). The uncertainty in 
A^87o(H2) was calculated by propagating the uncertainty in 7^"^ , 
and thus reflects only the Icr rms noise in the 870 //m map. 

The MSX 8 i^m optical thicknesses (Table H Col. (3)) were 
used to estimate the peak H2 column densities by applying the 
formula 



rdust 
'870 



fi870(7'd);UH2WHf870^d ' 



(8) 



A^8(H2) 



peak 



(9) 



where ct) is the dust extinction cross-section per H2 molecule. 
According to the OH94 dust model we have used (see Sect. 5.1), 
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Table 5. Results of the SED fits. 



Mo( -^^bol 7'cold T'hot M|o(/L(,ol 

Source [Mq] [10^ Lq] [K] [K] M,„m/M,„, L,„uIUo\ [Mq/Lq] 

IRAS 13037-6112 60 1.5 22 243 ~ 1 0.90 0.04 

IRAS 13039-6108 89 2.1 22 253 ~ 1 0.87 0.04 



Table 6. Linear radii, masses, and H2 column and volume-averaged number densities of the submm clumps within the IRDC 
G304.74. 









A?870(H2) 








<n(H2)> 


Source 


[pc] 


[Mq] 


[10^2 


cm 2] 


[W- 


cm--] 




[10-* cm-3] 


SMM r 


0.42 


107 


2.06 


±0.13 


0.74 


±0.30 


2.8 


± 1.1 


0.7 


SMM2[ 


0.33 


53 


1.35 


± 0.14 


0.54 


±0.27 


2.5 


± 1.3 


0.7 


SMM 33 


0.49 


182 


2.24 


±0.13 


0. 


14rf 


16 


±0.9 


0.7 


SMM4[ 


0.44 


219 


2.96 


±0.13 


0.36 


±0.24 


8.2 


±5.5 


1.2 


IRAS 13037~6112 


0.30 


48 


1.17 


±0.07 










0.8 


SMM 5 


0.30 


37 


0.85 


±0.13 


1.79 


±0.57 


0.5 


±0.2 


0.6 


SMM 6 


0.35 


96 


2.06 


±0.13 










1.0 


SMM 7 


0.35 


48 


1.17 


±0.14 


1.74 


±0.56 


0.7 


±0.2 


0.5 


IRAS 13039-6108 


0.47 


78 


1.18 


±0.08 










0.3 


SMM 8 


0.34 


59 


1.39 


±0.13 


2.23 


±0.76 


0.6 


±0.2 


0.7 


SMM 9 


0.37 


85 


1.70 


± 0.13 


3.13 


± 1.40 


0.5 


±0.2 


0.8 


IRAS 13042-6105 


0.31 


43 


0.94 


± 0.13 










0.7 



" By assuming = 30 K (see Sect. 4.3 and Tablegll, M,o„, = 39 Mq, WgyoCHj) = 0.75 ± 0.05 x 10" cm^^ A'gvo/A's = 1-0 ± 0.1, and 
<n(H2)) = 0.2 X 10* cm-\ 

* By assuming = 30 K, M,o„, = 20 Mq, A'87o(H2) = 0.49 ± 0.05 x 10^^ cm-^ A'gvo/A'g = 0.9 ± 0.1, and (^(Hz)) = 0.3 x 10"* cm-^. 
' By assuming Ti = 30 K, M^on, = 66 Mq, A'87o(H2) = 0.82 ± 0.05 x 10" cm^^ Ngio/Ns = 5.9 ± 0.4, and (^(Hz)) = 0.3 x 10"* cm^^. 
No error given because it is greater than the actual value. 

By assuming = 30 K, M^on, = 80 Mq, A'87o(H2) = 1.08 ± 0.05 x 10" cm-^ Ng^a/Ns = 3.0 ± 2.0, and (nCHj)) = 0.4 x 10* cm-^. 



the value cr^ - K^ifin^niuRd ~ 7.0 x 10 cm^ at 8.8 /urn. The 
uncertainty in A^8(H2) was calculated by propagating the uncer- 
tainty in Tg^j^'^ and do not include the systematic error resulting 
from the estimate of cta- The obtained H2 column densities, and 
the Nm)/Ns ratios, are given in Cols. (4)-(6) of Table|6] 

5.4. Extinction estimates witti 2MASS 

We examined if extinctions estimated from the JHR^^ photom- 
etry of 2MASS stars lying in the region can be used to cali- 
brate the H2 column densities derived abov^ A rarefaction of 
2MASS point sources can be discerned in the cloud region and 
its immediate vicinity, and it is not possible to derive an extinc- 
tion map of such detail that a comparison between the LABOCA 
map would be meaningful. Altogether 14 2MASS point sources 
with good photometric quality in all three bands lie within the 
LABOCA contour 0.1 Jy beam '. Three of them show J - H 
and H — K^ colours characteristic of giant stars reddened by sub- 
stantial columns of interstellar dust. This judgement is based on 
their locations near the standard reddening line on the J - H vs. 
H - Ks plot, and on the fact that their J - excesses, E j^k, are 
larger than 0.6 (see below). The three stars are the 2MASS PSC 
objects 13064180-6128529, 13064266-6128213, and 13070908- 
6124303 (see TableQand Fig.|6l). The first two lie near the clump 
SMM 4 and are likely to represent background K or M giants. 
The third is found very close (~ 10") to IRAS 13039-6108 and 
is somewhat too bright for a class III giant in view of the distance 
and extinction. Because its NIR colours (J - H - 2.56 + 0.03, 
H — - 1.31 + 0.03) are consistent with an YSO (young stel- 



' The 2MASS All-Sky Point Source Catalog (PSC) used here have lar object) candidate with NIR excess (Matsuya nagi et al. 2006) l, 
been made available at h ttp : //irsa . ipac ■ caltech . edu/i the source is possibly a NIR counterpart of IRAS 13039-6108. 




4 2 -2 -4 

Aa (arcmin) 



Fig. 6. LABOCA dust continuum map (contours) with 2MASS 
point sources (plus signs). The heavily reddened 2MASS stars 
within the cloud boundaries are marked with diamonds. The con- 
tour levels go from 0.1 to 0.6 by 0.1 Jy beam"'. 
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J 


H 




J-K, 




-'V2MASs(H2!l^ 

[10^2 cm-2] 


-'V87o(H2) 




2MASS designation 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[10^2 cm.--] 


A'87o/A'2MASS 


13064180-6128529 


14.15 ±0.04 


12.79 ±0.05 


12.04 ± 0.04 


2.11 ±0.06 


1.0- 1.5 


0.53 - 0.80 


2.3 


2.9-4.3 


13064266-6128213 


15.98 ±0.08 


14.42 ± 0.04 


13.53 ± 0.05 


2.44 ±0.10 


1.3-1.8 


0.71-0.98 


1.4 


1.4-2 


13070908-6124303 


14.26 ±0.02 


1 1.70 ± 0.02 


10.39 ± 0.02 


3.87 ± 0.03 


2.7-3.2 


1.5 - 1.7 


0.62 


~ 0.4 



" Assuming (J - A')o = 0.63 - 1.13 corresponding to spectral classes KO III-M3 III. 
* Assuming N{H2)/Ej-k = 5.4 x 10^' cm^^ mag"'. 



The H2 column density ranges given in Col. (7) of 
Table |7] are derived using intrinsic colours for KO III- 
1VI3 III giants (Bessell & Brett 1988), and the relationship 
A^(H2)/£'j-^ = 5.4 X 10^'cm '^mag'' (Bohli netal. 19781 
IJVIathis 19901 [Harjunpaa & Mattila 1996, 1. The colour excesses, 
E j-K, measure the total amount of dust in front of the stars, so 
they are affected by dust in the foreground, and possibly also 
in the background of the cloud. Therefore, the A^(H2) values are 
upper limits for the contribution of the cloud itself. A nearby 
line of sight {I = 304?75, b = \°.25) has been included in the 
investigation of 3-D distribution of the extinction by Marshall 
et al. (2006). Their results suggest a foreground extinction of 
Ak ~ 0.4 (Ej^K ~ 0.6; Marshall et al. 2006..Mathis 1990J up to 
the cloud's distance 2.4 kpc. 

The LABOCA 870 fim intensities in the directions 
of 13064180-6128529, 13064266-6128213, and 13070908- 
6124303 are 0.49, 0.31, and 0.13 + 0.03 Jy beam"', respectively. 
The conversion to the H2 column densities using Eq. ^ with 
Td = 15 K gives 2.3 x lO^^, 1.4 x 10^2, and 0.62 x 10^^ cm^^. 
The first two values are not consistent with the upper limits de- 
rived from Ej^K (see Tabled Col. (9)). In the case of 13064180- 
6128529, the H2 column density derived from LABOCA 870 fim 
data can be brought down to conform with 2MASS estimate by 
increasing the dust temperature to = 30 K, whereas for the 
location of 13064266-6128213 an increase to I'd = 20 K would 
be sufficient. The latter value is in good agreement with the bolo- 
metric temperature of IRAS 13037-6112 (Tcoid - 22 K, see 
Table |5j Col. (4)) which lies (in the plane of the sky) quite close 
to 13064266-6128213. The 8 fim intensities in the directions of 
13064180-6128529 and 13064266-6128213 are 1.03 x 10"'' and 
1.22 X 10"^ W m"2 sr"', respectively. These coiTespond to 8 
fim optical thicknesses of 0.41 + 0.19 and 0.02 (here the asso- 
ciated eiTor is larger than the value), respectively (see Eq. (|2]|). 
Using Eq. the corresponding H2 column densities become 
0.59 ± 0.27 X 10^2 and 0.03 x 10^^ cm^^, respectively. The first 
value is consistent with the 2MASS estimate, whereas the latter 
value is much lower. 

The above results clearly show that it is important to know 
the dust temperature in order to accurately determine the H2 col- 
umn density from the submm dust continuum emission (Eq. (O). 
On the other hand, the MIR absorption and 2MASS extinction 
methods require several uncertain assumptions, such as the dust 
model and the corresponding dust extinction cross-section, con- 
tribution of the foreground emission, and the relation between 
A^(H2) and Ej^k- Moreover, as the 2MASS extinction could be 
estimated only along a few lines of sight, the present statistics is 
very poor Thus, the data presented here do not offer a firm con- 
clusion about the most reliable method to determine the H2 col- 
umn densities. We note that the empirical relationship between 
A^H and A j derived by Vuong et al. (2003) in p Oph, used recently 
by Marshall et al. (2009), implies N{H2)/Ej^k ^ 4.7 x 10^' 



cm 2 mag ^ This would make the H2 column densities about 



13% smaller than those estimated above. 



6. Discussion 

6.1. Estimating tiie background and foreground MIR 
emission 

The average MIR emission around the cloud, /mir = 1.23+0.02x 
10~* W m"2 sr was estimated in Sect. 4 by using median fil- 
tering. Vasyunina et al. (2009) estimated /mir in the vicinity of 
several IRDCs directly from the 8 fim images. They pointed out 
that it is difficult to control the process of median filtering if 
large (> 10') filters are needed. On the other hand, if the filter is 
too small (comparable to the cloud size), the value of /mir will 
be underestimated because it will be affected by the cloud itself 
(Butler & Tan 2009). Because the largest diameter of G304.74 
is ~ 13', we tried the direct method by choosing manually MIR 
emission patches (free of strong MIR emission sources) in the 
close vicinity of the IRDC. We used three different patches with 
angular sizes of 6'.'7 x 3'.'3, 4'.'8 x 2'.'3, and 5'.'7 x 1'.'8. The mean 
and standard deviation of the 8 fim intensity within these regions 
were 1.17 + 0.18 x lO"'', 1.21 +0.18x 10"^ and 1.25 + 0.21 x 10"^ 
W m"2 sr respectively. The mean and standard deviation of 
these three are 1.21 ± 0.11 x 10 W m"^ sr"'. This result is in 
excellent agreement with the value obtained from median filter- 
ing. 

Based on the observed anti-correlation between the 850 fim 
and 8 fim intensities (cf. Fig. HJ, Johnstone et al. (2003) and 
Ormel et al. (2005) found that /foie - he. in the case of IRDC 
Gil. 11-0.12 and the W51 IRDC, respectively. Peretto & Fuller 
(2009) constrained /fore by the requirement that both MIR ab- 
sorption and 1.2 mm emission should give the same A^(H2). In 
this method, it is assumed that the true 8 fim opacity can be cal- 
culated from the millimetre emission (see Eq. (3) of Peretto & 
Fuller (2009)). Then, the value of /fore can be calculated in terms 
of /mir (see our Eq. (|2]i). Peretto & Fuller (2009) also showed 
that on the average the background emission is approximately 
equal to the foreground emission (their Eq. (5)). The drawback 
in this approach is the uncertainty in the ratio of the mass ab- 
sorption coefficients at MIR and (sub)mm wavelengths (R/^ in 
Eq. (3) of Peretto & Fuller (2009)). In the present study, we esti- 
mated the contributions of background and foreground emission 
using the foreground intensity ratio (Sect. 4.2). We also ended 
up with the result that /fore - hg- Butler & Tan (2009) and 
RBG09 estimated that /fore - (0.1 - 0.5) x /bg (/fore-method) 
and /fore = (2 - 5) X /bg (850 ^(m-8 fim anti-correlation) for their 
samples of IRDCs, respectively. 
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6.2. Comparison of H2 column densities determined from 
dust continuum emission and extinction data 

The H2 column densities derived from submm emission and 
MIR absorption are mostly in good agreement (within a factor 
of ~ 2; Table|6] Cols. (4) and (5)). This suggests that the dust pa- 
rameters used in Eqs. ([8]l and ^ are reasonable. The dust tem- 
perature, Td, is likely to show spatial variations, and this causes 
uncertainties to column density estimates based on submm emis- 
sion (Eq. (I8]l), but does not affect the 8 /^m absorption method 
(Eq. The greatest differences between the two A^(H2) val- 
ues are found toward positions with the highest A^87o(H2) values. 
These are SMM 1, 3, and 4, where the A''87()(H2)/A''8(H2) ratios 
are 2.8 ± 1.1, 16 ± 0.9, and 8.2 ± 5.5, respectively (Col. (6) of 
Table |6]l. Some of these discrepancies could (in part) be reme- 
died, e.g., by increasing the dust temperature (cf. footnote in 
Table |6]l. It seems possible that the dust temperature is higher 
than 15 K in the four clumps (SMM 1^) near southwestern tip 
of the cloud. On the other hand, Vasyunina et al. (2009) showed 
that the extinction method becomes unreliable at very high col- 
umn densities, but this should happen only at A^(H2) > lO^^* 
cm"^. 

However, there are also uncertainties in deriving A^8(H2), in 
particular related to the contamination by the foreground emis- 
sion (here done by using the foreground intensity ratio, see 
Sect. 4), and the dust extinction cross-section (reliable only 
within a factor of ~ 2, e.g., Ragan et al. 2006). Also, the value of 
ka, needed in the calculation of A^87o(H2) (Eq. ([8]l) has an uncer- 
tainty similar to that of cr^. The dust-to-gas mass ratio could also 
differ from the adopted value 1/100 as mentioned in Sect. 5.1. 
Thus, the direct comparison of A^(H2) derived from dust contin- 
uum and extinction data should be taken with caution. We note 
that Parsons et al. (2009) found for their sample of IRDCs that 
Ns{U2) (MSX) and A^85o(H2) (SCUBA) agree within an order 
of magnitude. The moderate correspondence may be partly ex- 
plained by the omission of the foreground emission. 

6.3. Nature of submm clumps 

The MSX 8 fim emission associated with molecular clumps 
suggests the presence of protostars. By combining the submm 
LABOCA and MIR MSX data one can distinguish between can- 
didate starless and star-forming clumps. Of the twelve clumps in 
G304.74, four are likely to be associated with newly born stars. 
In addition to the three IRAS (and MSX) sources, the clump 
SMM 6 is associated with two MSX 8 fim point sources. The 
remaining eight MIR dark clumps are either starless or con- 
tain low-luminosity (low-mass) YSOs falUng under the detec- 
tion limit of the MSX (cf. Parso ns et al. 2009b . For instance, the 
Spitzer/GLIMPSE survey with better sensitivity and resolution 
compared to that of MSX, revealed that some IRDC clumps 
previously thought to be starless do contain MIR sources (e.g., 
[Chambers et a l. 2009). Moreover, Chambers et al. (2009) de- 
tected H2O masers towards some MIR dark clumps, which is 
a clear indication of star formation activity in them. On the other 
hand, starless clumps may be either prestellar or just unbound 
structures that will eventually disperse. 

The extremely short lifetime of starless IRDC clumps (10^ - 
10"* yr, i.e., ~half of the time spent in the protostellar phase) 
derived by Parsons et al. (2009) does not support the idea that 
all of the eight dark clumps in G304.74 could be starless. If 
gravitationally bound, these clumps with masses in the range 
~ 40 - 200 M p (Table E Col. (3)) are capable of forming 
high-mass stars ( Thompson et al. 2006 1, and some of them could 



represent/contain HMSCs (cf. [Chambers et al. 20091) . The peak 
column densities of the clumps do not reach the minimum col- 
umn density threshold of 1 g cm"^, i.e., A^(H2) ~ 2.2 x lO^-' 
cm"^, proposed by Krumholz & McKee (2008) for the forma- 
tion of high-mass stars. However, it is possible that the clumps 
host cores where such a high A^(H2) values can be reached (e.g., 
IHennemann et al. 20091) . 

The presence of high-luminosity (~ 1.5 - 2 x 10^ 
L0) IRAS sources indicates that intermediate to high-mass 
star formation is going on in other parts of the cloud. 
In fact these luminosities suggest intermediate-mass stars as 
high-mass protostellar objects (HMPOs) are often found in 
the range ~ l O "^-^ - 10^-^ Lq (e.g. Sridharan et al. 20021 
IFazal et al. 20081 IGrave & Kumar 2009I I. On the other hand, 
the envelope masses of these sources (~ 60 - 90 M0, see 
Col. (2) of Table |5j are sufficiently large for high-mass star 
formation dBeuther & Steinacker 2007 ). For comparison, major- 
ity of the low-mass protostars in nearby (d < 250 pc) molec- 
ular clouds have luminosites Lboi ^ a few x Lq, the high- 
est observed values being ~ 70 Lq in only a few cases (e.g., 
Evans et al. 2009). IRAS 13037-6112 and IRAS 13039-6108 
have similar SED parameters (see Fig.|5]and Table|5j, suggesting 
that they probably represent the same evolutionary stage. Based 
on its IRAS [25-12] and [60-12] colours (0.61 and 1.55, respec- 
tively), IRAS 13037-61 12 belongs to the so-called High sources, 
and could potentially be associated with ultra-compact (UC) 
HIT region (e.g., [ Wood & Churchwe ll 19891 IPallaetal. 19911 
IMolinari et al. 1996) . However, the M/L ratio of IRAS 13037- 
6112 (0.04) is more typical of sources younger than UC HII re- 
gions jSridharan et al. 2 002 ). IRAS 13039-6108 has a [25 - 12] 
colour index of 0.27, settling it to the so-called Low sources. 
Moreover, Fontani et al. (2005) found that the C'^0 and CS 
linewidths in IRAS 13039-6108 are significantly smaller than 
those typically observed in massive clumps associated with 
UC HII regions (e.g., ICesaroni et al. 19911 IHofner et al.lOOOl l. 
This further supports the idea that both IRAS 13037-6112 and 
IRAS 13039-6108 represent the same evolutionary stage (ear- 
lier than UC HII). Using the assumption of an isothermal sphere 
with Tkin = 22 K and a density profile of the form n(r) oc 
r"'-^, which is typical of high-mass star-forming clumps (e.g., 
IBeuther et al. 2002al l. the C'^0(2 - 1) linewidth of 0.93 km s"' 
observed by Fontani et al. (2005) imphes a virial mass of ~ 99 
M0 for IRAS 13039-6108 (M^ont/M^i, ^ 0.8, MsEo/Myt ^ 0.9, 
i.e., the virial parameter defined by Bertoldi & McKee (1992) 
is gyir = M^iJM ^ 1.1 - 1.3) (see, e.g., Eqs. (1) and (2) 
of I Chen et al. 2008l l. Thus the clump is near virial equilibrium. 
IRAS 13042-6105 in the northern part of the cloud is probably 
in an earlier stage of evolution than the other two IRAS sources 
because it is not as bright at FIR wavelengths (see Table [3]i. 

6.4. Clump mass distribution 

The mass distribution of clumps/cores is important parameter 
concerning the cloud fragmentation mechanism. Our sample of 
clumps is, however, so small (12 in total, 8 MIR dark) that 
it is not reasonable to study their mass distribution directly. 
Therefore, we only compared it with the mass distributions de- 
rived for other, larger IRDC clump samples, using the results 
of Sridharan et al. (2005), RJS06, Vasyunina et al. (2009), and 
RBG09. 

Figure I2] presents the observed cumulative mass functions, 
which include clumps of mass less than M, i.e., N{M) - N{m < 
M)/Ntoi, for G304.74 and for a sample of cold clumps from 
RJS06 and RBG09. We note that the IRDC clump mass func- 



Miettinen & Harju: LABOCA mapping of the IRDC G304.74 



11 



tions in RJS06 and in our work are constructed by removing the 
MIR bright clumps from the samples (in the case of G304.74, 
this means the three IRAS sources and SMM 6). From the sam- 
ple of Sridharan et al. (2005), we excluded the high temperature 
clump No. 9, because its temperature (32.7 K) was much higher 
than the rest of the sample. From the the sample of RBG09, we 
removed the clumps associated with YSOs, and the clumps pos- 
sibly contaminated by foreground (or background) stars. Thus, 
the mass functions include only those clumps that initially have 
all their mass available for star formation. 

The previous studies taken into this comparison used slightly 
different assumptions about the dust temperature and opacity. 
We scaled all the clump masses included in the cumulative mass 
functions to correspond to a uniform dust temperature of Td = 15 
K, and an opacity that is consistent with our k^jq = 0. 17 m^ kg"' . 
RBG09 derived clump masses from the total hydrogen column 
densities, A^h, as estimated from the 8 //m optical thicknesses, rg. 
For the dust extinction cross-section per H nucleus at 8 //m, erg, 
they used the value 2.3 x 10"^^ cm^, based on the Weingartner 
& Draine (2001) dust model. According to the model we have 
used ( jOssenkopf & Henning 1994) l, the corresponding number 
is 3.5 X 10"^^ cm^ per H nucleus (see Sect. 5.3). In their mass 
formula, RBG09 used a factor of 1 . 16 as the ratio of the total gas 
mass (including He) and the hydrogen mass. In our calculations 
this ratio has been 1.4 (Sect. 5.3). These differences have been 
accounted for by scaling the masses from RBG09 by 0.8 in the 
comparison with our results. 
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Table 8. K-S test results between different IRDC clump mass 
distributions. 
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Fig. 7. Normalised cumulative mass functions, N(M), for the 
MIR dark clumps in the IRDC G304.74 (solid line) and in sev- 
eral other IRDCs studied by RJS06 (dotted line), and by RBG09 
(dashed line). Clump masses from RJS06 and RBG09 were 
scaled down to match our assumptions about and /c^ (see main 
text) 



To determine whether our clump masses and those from 
other studies are derived from the same clump mass distribution, 
we carried out the two-sample Kolmogorov-Smirnov (K-S) test. 
For this test, the mass scales were matched, i.e., the comparison 
was done within the range of common mass interval. The K-S 
test results are shown in Table |8] The columns of this table are: 
(1) the survey used in the comparison; (2) number of clumps 
included; (3) the maximum vertical difference between the cu- 
mulative mass distributions (Z)max); (4) the probability for the 





Number of 






Study 


cold clumps 


^max 


Prob. 


Ragan et al. (2009, RBG09J;;, 


345 


0.153 


0.998 


Rathbome et al. (2006, RKObf 


140 


0.272 


0.739 


Sridharan et al. (2005) 


28 


0.359 


0.534 


Vasyunina et al. (2009) 


17 


0.449 


0.408 



" The probability that the RBG09 and RJS06 clump mass sam- 
ples represent the subsamples of the same parent distribution is 100% 
(£»„,„ = 0.001). 



null hypothesis that the two functions are drawn form the same 
parent distribution (the significance level of the K-S statistic). 

The highest number statistic is provided by the study of 
RBG09, and when compared with this study, the K-S test yields 
a probability of 99.8% that the clump mass distributions in 
G304.74 and the other IRDCs are drawn from the same parent 
distribution. This probability drops significantly when smaller 
samples are used in the test (Table [8] Col. (4)). Moreover, ac- 
cording to the K-S test, there is a 100% probability (D^ax = 
0.001) that the RBG09 and RJS06 clump masses represent the 
subsamples of the same underlying parent distribution. For this 
test, clump masses from RBG09 were also multiplied by 0.8 to 
compare with the RJS06 masses which were based on the value 

f 1.2mm = 0.1 m^ kg"'. 

For the clump masses between ~ 30 and 3000 Mq, RBG09 
derived an IRDC clump mass spectrum with a slope of a = 
1.76 ± 0.05, which is consistent with the mass functions derived 
for high-mass star-forming regions, and also resembles the mass 
function of Galactic stellar clusters (see RBG09 and references 
therein). Note, however, that RBG09 merged all their clumps 
into a single mass function, whereas we have used only starless 
clumps in deriving the cumulative mass functions. RJS06 found 
Salpeter 1955[ l mass function for their 



a Salpeter-like (a = 2.35; 

cold IRDC clumps (a = 2.1 ± O.TatM > 100 M©). The IRDC 
clump mass spectra derived by RJS06 and RBG09 are compa- 
rable to those predicted by turbulent fragmentation models (see 
Sect. 6.7 and references therein). 

6.5. Clump spatial distribution 

In addition to the clump mass distribution, it is also useful to de- 
termine how clumps are spatially distributed within the cloud in 
order to unveil the presence of a possible preferred length-scale 
of fragmentation (e.g.. lMufloz et al. 2007I I. For this purpose, we 
determined the clump-separation distributions and the number 
distributions of the projected separation distance between near- 
est neighbour^ in G304.74, and in nine other IRDCs for com- 
parison. We chose those IRDCs from the sample of RJS06 
which contain the largest number of clumps, i.e., MSXDC 
G023.60h-00.00, G024.33-h00.11, G028.37-h00.07, G028.53- 
00.25, G03 1.97+00.07, G033.69-00.01, G034.434-00.24, and 
G035. 39-00. 33. Moreover, we determined the spatial distribu- 
tion of YSOs in the IRDC MSXDC G048.65-00.29 studied by 
van der Wiel & Shipman (2008). For these analyses, we selected 
only those sources that are clearly associated (in the plane of the 
sky) with their parental dark cloud (e.g., clumps that lie within 



If the number of clumps in the cloud is A'ci, then the number of 
clump-separations is A'ci (A'ci - l)/2, whereas the number of distances 
between an individual clump and its nearest neighbour is equal to A'ci. 
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the dark filaments). Thus, we excluded the millimetre clumps 
MM 1, 3, and 5 from G023.60+00.00, MM 2, 5, and 7 from 
G024.33+00.11,MM3,5,7,8, 12, 13, 18 from G028. 37+00.07, 
MM 1 and 2 from G035.39-00.33, and YSOs S4, 10, 17, 18, 19, 
and 20 from G048. 65-00. 29. In addition to the observed spatial 
distributions, we also determined the distributions expected from 
random positions of the same number of objects as the observed 
samples have. The areas over which the objects were randomly 
distributed were chosen so that they approximate the observed 
dark cloud areas; the IRDC areas were estimated by rectangles 
which just cover the observed dark clouds. When needed, these 
rectangles were rotated with respect to the {a, 5)-coordinate sys- 
tem. The random distributions were generated a hundred times 
per cloud and the resulting averaged histograms were used in 
comparisons with observed spatial distributions. 

Figure [8] (top) shows the observed clump-separation distri- 
bution in G304.74, and the distribution expected for the same 
number of randomly positioned clumps over minimum rectan- 
gular area which encloses the dark cloud (~ 29.8 arcmin^). 
The mean and its standard deviation, and median of the clump 
separations in G304.74 are log«r)obs/AU) = 5.690 + 0.041 
(4.90!JJ;4^ X 10^ AU) and log(robs/AU) = 5.759 (5.74 x 10^^ 
AU), respectively. These values are similar to those of ran- 
domly positioned clumps, for which the mean and median are 
log«r)ra„/AU) = 5.674+0.061 andlog(rra„/AU) = 5.738+0.084 
(see Table |9]l. The latter two values and their +-errors quoted 
represent the average values and their standard deviations de- 
rived from the 100 generated random distributions. According 
to the K-S test, the probability that the observed distribution and 
the generated random distribution represent the same underlying 
distribution is 100%. Statistics of the clump-separation distri- 
butions in other IRDCs studied by RJS06 and van der Wiel & 
Shipman (2008) are listed in Table |9] The columns of this table 
are the following: (1) IRDC name; (2) number of clumps used 
in the analysis (see above); (3) distance; (4) area used to cre- 
ate the random distributions (see above); (5) and (6) mean and 
median of the observed clump-separation distribution ({r)obs and 
'^obs); (7) and (8) mean and median of the corresponding random 
distribution ((r)i-an and rran); (9) and (10) ratios between the ob- 
served and random mean and median separations (quoted errors 
are propagated from the standard deviations of (r) and r); (11) 
probability given by the K-S test that the observed and random 
distributions are drawn from the same underlying distribution. 
The observed clump separations are mostly similar to those ex- 
pected from random distributions. This is evident from the ra- 
tios (r)obs/('")ran and Tobs/'^ran which are close to unity, and from 
the K-S probabilities which are high (~ 71 - 100%) except for 
three cases (G035.39, G028.37, G024.33; see Cols. (9)-(ll) of 
Tableig. 

Figure |8] (bottom) compares the observed nearest-neighbour 
distribution in G304.74 with the distribution for randomly posi- 
tioned clumps. The mean and median of the nearest-neighbour 
distribution in G304.74 are log«r)obs/AU) = 5.083 + 0.058 
(1.21+;j [^ X 10-** AU) and log(robs/AU) = 5.136 (1.37 x 10^ 
AU), respectively. Again, these values are comparable to those 
of randomly positioned clumps, for which the mean and me- 
dian are log«r>ran/AU) = 5.030 + 0.111 and log(r,-an/AU) = 
5.047 + 0.1 10, respectively (see Table[TOll. According to the K-S 
test, there is about 90% probability that the observed and random 
nearest-neighbour distributions are samples of the same under- 
lying distribution. We note that the minimum observable sepa- 
ration corresponds to the beam size, i.e, 18'.'6 or ~ 4.46 x lO'* 
AU (log(r/AU) = 4.649) at 2.4 kpc. Statistics of the nearest- 
neighbour distributions in other IRDCs are given in Table [TOl 



Table [TOl have the same meaning as in Table |9] except now for 
nearest neighbour separations. The observed nearest-neighbour 
distances are similar to those expected from random distribu- 
tions. This is evident from the ratios (r)obs/('")ran and robs/'^ian 
which are (within the erros) about 1, and by the high K-S proba- 
bilities (~ 59-100%) in every other case except G03 1.97, where 
this probability is still ~ 37% (see Cols. (6)-(8) of Table [TOll. 

In summary, the average projected separations between 
clumps in the studied IRDCs range from about 2.6 x 10^ AU 
to 1.2 X 10^ AU (i.e., the minimum and maximum lie within a 
factor of five), and the average projected distances between the 
nearest neighbours range from 6.0 x 10'* AU to 3.5 x 10^ AU 
(i.e., the changes are within a factor of six). For most clouds, 
the distributions of projected separations and distances between 
the nearest neighbours can be mimicked by clumps placed ran- 
domly into the same projected area as occupied by the cloud. 
Assuming that the vectors connecting clump pairs are randomly 
oriented, the average projection factor is (sin = n/4, where 
dij is the angle between the line of sight and the vector pointing 
from clump / to clump j. Correcting for this projection effect, 
the grand averages of the separations and distances between the 
nearest neighbours are about 6.5 x 10^ AU (3.1 pc) and 2.2 x 10^ 
AU (1.1 pc), respectively. 

6.6. Fragmentation of IRDCs 

One plausible scenario for the origin of filamentary clouds 
is that they are formed in shocks occurring in converging 
flows driven by large-scale turbulence (e.g.. lKlessen et al. 20001 
Pad oan et al. 2 001). The chaotic process can give rise to ran- 
domly positioned density peaks within filaments, and these can 
become centres of gravitational collapse. On the other hand, sup- 
posing that compression leads to an equilibrium structure, a fil- 
ament can fragment through the Jeans instability. 

The critical wavelength, A^, of perturbations leading to grav- 
itational instability depends on both the gas kinetic tempeture, 
Tkin, and the density, p: ~ c^/ ^jGp, or, in terms of the surface 
density, E: /Ic ~ c^/CL, where Cs is the sound speed, and G is the 
gravitational constant (e.g., Larson 1985IIHartmann 2002l l. 

The determination of the 'Jeans length' is not quite straight- 
forward in a study based on dust emission because the cloud 
mass and therefore also the average density and surface density 
depend on the assumed dust temperature, Tj. Furthermore, in the 
case of G304.74, there is no independent estimate of Tkin, but it 
is assumed to be equal to T^. 

The assumption = 15 K yields a total mass of ~ 1000 
Mq and an average surface density of 0.05 g cm"^ for G304.74 
(within the LABOCA contour 0.1 Jy beam '). In these cir- 
cumstances, the critical wavelength in an isothermal equilib- 
rium filament is = 0.19 pc or 40000 AU ( IHartmann 20021 
Larson 1985), and the corresponding mass is Mc ~ 5 Mq. The 
comparison between dust emission at 870 yum, 8 //m absorption, 
and the visual extinction from 2MASS suggests an elevated tem- 
perature in the southern part of the cloud (Table|4] Col. (4)). By 
assuming = 30 K, one would obtain a total cloud mass of 
~ 400 Mq, and an average surface density of E = 0.02 gcm"^. 
These values of temperature and surface density imply A^ = 1.1 
pc or 2.2 X 10^ AU, and ~ 50 Mq. The critical wave- 
length, Ac, is not expected to determine a uniform length scale 
of fragmentation. According to the analysis of Stoddlkiewicz 
(1963; see also Curry 2000 and references therein), the fastest 
growing perturbations have a length scale of roughly twice A^. 
Nevertheless, the fragmentation of a homogenous cloud is likely 
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to result in a preferred length scale and a quasi-periodic struc- 
ture, as opposed to a random distribution of clumps. 

The projected distances between the nearest neighbours in 
G304.74 lie in the range 4 x 10"^ - 2.4 x 10^ AU. They are 
comparable to the characterics length scales indicated above. 
However, the clump masses in the cool northern part (~ 40 - 90 
Mq) are about ten times larger than Jeans masses at 15 K and 
the surface density implied by this temperature. Either the Jeans 
instability has occurred when the cloud has been warmer, i.e. 
~ 30 K, and the clumps have cooled during the contraction, or 
the dense filament is a result of strong compression by external 
forces. We note that the filament is thinner at the northeastern 
end (cross-sectional diameter ^ ~ 0.45 pc) than in the south- 
west ((p ~ 1 .4 pc). Moreover, the filament can have accumulated 
mass through gravitational inflow from the surrounding cloud 
dHeitsch etal. 2009l l. 

The present average mass line density in G304.74, ~ 100 
M0 pc"' exceeds the critical value for a non-magnetic, self- 
gravitating isothermal cylinder in equilibrium (~ 25 M0 
pc"' at 15 K, the critical line density directly proportional 
to Tkin; lOstriker 19641 [Curry 2000| . For comparison, in sev- 
eral filamentary IRDCs studied by Rathborne et al. (2006; 
G025.04, G028.53, G028.53, G031.97, G033.69, G034.43, and 
G035.39), the Une densities are in flie range ~ 70 - 800 M© 
pc"'. Supercritical line densities seem to be a common fea- 
ture in filamentary IRDCs, and they are likely to fragment 
into smaller cores. Our observations cannot resolve the pos- 
sible fragments in G304.74, but the process is manisfest in 
presence of (intermediate- to high-mass) protostars, i.e., the 
two IRAS sources. On the other hand, the recent interfero- 
metric studies of G28. 34 (Zha ng etal.2009| and IRDC 19175 
( Beuther & Henning 2009) have provided direct evidence for 
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sub-fragmentation of IRDC clumps. The two studies offer, how- 
ever, contradictory views of the nature of the underlying insta- 
bility (turbulent fragmentation vs. Jeans instability). 

6. 7. The origin of IRDCs and substructures witliin them 

IRDCs are the densest parts of molecular cloud complexes. 
Several observational facts and modelling results suggest 
that turbulence has an important role in the formation of 
IRDCs and their fragmentation into clumps. 1) On large 
scales, molecular clouds have highly supersonic linewidths. 
In general, these are considered to imply turbulent mo- 
tions (e.g., [Elmegreen & Scalo 2004[ iMcKee & Ostriker 2007l i. 
However, large linewidths in dense clouds can also be ex- 
plained by collapsing motions towards local gravitational cen- 
tres dHeitsch et al. 2009l l. 2) The filamentary shapes of IRDCs, 
and molecular clouds in general, are consistent with cloud 
morphologies predicted by numerical models of supersonic 
turbulence driven on large scales (e.g., IKlessen 20001 2001; 
[Jappsen et al. 2005)1. In this model, a dense filament can form 
where converging flows meet. Also the mass surface density 
distribution observed in IRDCs have been found to correspond 
to expectations from numerical simulations of turbulent clouds 
dButler & Tan 2009] l. 3) The high-mass end of the mass dis- 
tiribution of a large sample of IRDCs studied by Marshall et 
al. (2009) follows a power-law {AN I AM cc M"'^^) which can 
be reproduced by density fluctuations induced by turbulence. 
4) The mass spectrum of the high-mass clumps within the ex- 
tensive IRDC sample of RJS06 can be fitted with a power- 
law {AN/AM oc M"^ '). This power-law agrees with fragmen- 
tation due to supersonic turbulence in self-gravitating clouds 
dKlessen 20010 . 5) Finally, the spatial distributions of clumps 
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Fig. 8. Top: Observed clump-separation distribution (solid line) 
compared with the expected distribution for random distribution 
of the same number of clumps as the observed sample over rect- 
angular area which approximate the area of the IRDC (dashed 
line). Bottom: Observed nearest-neighbour distribution (solid 
line) compared with the expected distribution for random dis- 
ti'ibution (dashed fine). 



within IRDCs (cf. Sect. 6.5) show no clear deviation from a ran- 
dom distribution. This is what can be expected if fragmentation 
is driven by a stochastic process. 

The filamentary stuctures of IRDCs and the fact that star 
formation in them takes place in clusters suggest that tur- 
bulence is driven on large scales (e.g., Klessen 2001). The 
most likely driving agent for large-scale turbulence is pro- 
vided by supernova explosions, with an important contribu- 
tion of density fluctuations caused by older remnants (e.g., 
Korpi et al. 1998; Joung et al. 2009 and references therein). 
The Galactic distribution of IRDCs peak in the 5-kpc 
molecular ring dSimon et al. 2006al [Marshall et al. 20091 1, where 
most of the Galactic supernova remnants are found (e.g., 
Jack son et al. 20081 1. 

The role of turbulence in the fragmention of IRDC clumps 
into dense cores is less evident. The observed spectral linewidths 
in clumps are broader than expected from thermal broadening 
(e.g., 'Ragan et al. 20061 ISakai et al. 2'008l IGibson et al.l009] l. 
but this does not inevitably imply that the gas is turbulent 
dHeitsch et al. 2009] ). As discussed in Sect. 6.6, high-resolution 
studies have revealed some regular stuctures within clumps 
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Table 9. Statistics of the clump-separation distributions in IRDCs. 
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Table 10. Statistics of the number distributions of the projected separation distance between nearest neighbours in IRDCs. 
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G024.33+00.11 


5.236 


+ 0.040 


5.230 


5.308 + 


0.117 


5.314 + 





148 


0.85 


+ 0.24 


0.82 + 


0.28 


0.792 


G023.60+00.00 


5.280 


+ 0.057 


5.341 


5.142 + 


0.148 


5.153 + 





172 


1.37 


+ 0.50 


1.54 + 


0.61 


0.780 



( [Zhang et al. 2009[ [Beuther & Henning 2009) 1 which point to- 
wards gravitational or fluid dynamical instabilities. In any case, 
self-gravity of shock-compressed clumps or filaments is inherent 
also in the models of turbulent fragmentation. 

7. Summary and conclusions 

We have mapped IRDC G304.74 in the 870 //m dust continuum 
emission with the APEX telescope. The submm dust continuum 
observations have been used together with 1.2 mm data from 
Beltran et al. (2006), and archival MSX and IRAS infrared data 
to derive the physical characteristics of the clumps within the 
cloud. Besides the dust continuum we used dust extinction data 
from MSX and 2MASS to derive the H2 column densities, and 
the mass distribution in the cloud. The results obtained via differ- 
ent methods are in reasonable agreement with each other taking 
into account the uncertain nature of some of the dust proper- 
ties and the relation between the H2 column density and extinc- 
tion. However, the agreement can be improved by assuming an 
elevated temperature in four clumps near the southwestern end 
of the cloud, and the possibility of a temperature gradient from 
about 15 K in the north and centre to about 20-30 K in the south 
cannot be ruled out. 

The filamentary cloud G304.74 contains 12 submm clumps. 
Star formation has already started in the cloud as three of the 
clumps are associated with both MSX and IRAS point sources. 
The SEDs of the two IRAS sources indicate bolometric lumi- 
nosities in the range ~ 1.5 - 2 x lO^* Lq. These are likely 
to be intermediate- or high-mass protostars. In addition, one 
of the clumps (SMM 6) is associatied with two MSX 8 yum 
point sources. The remaining eight submm clumps are MIR 



dark. The masses of these clumps (~ 40 - 200 M0) are suffi- 
ciently large to enable high-mass star formation, or alternatively, 
they can represent the cold precursors of proto-clusters. Thus, 
some of the candidate starless clumps in G304.74 could rep- 
resent/harbour the so-called high-mass starless cores (HMSCs, 
e.g., IBeuther et al. 2007] l. Further studies of these high-mass 
starless clump/core candidates are important in order to con- 
strain the initial conditions of high-mass star and star cluster 
formation. 

The clump masses in G304.74 were compared with the 
clump mass spectra from more extensive surveys of IRDCs. We 
found that IRDC clump masses from the present work, and those 
from Rathborne et al. (2006) and Ragan et al. (2009) probably 
represent subsamples of the same parent distribution. Also, av- 
erage distances between a clump and its nearest neighbour in 
different IRDCs are comparable (within a factor of ~ 6), sug- 
gesting that the fragmentation length-scale does not vary much 
from cloud to cloud. Moreover, in most IRDCs, clumps seem to 
be randomly distributed within the cloud area. These character- 
istics, and the fact the star formation in IRDCs predominantly 
occurs in the cluster mode, agree with models where fragmenta- 
tion is driven by large-scale turbulence. It is not clear, yet, how 
effectual turbulence is for the fragmentation of IRDC clumps 
into dense cores. High spatial resolution studies, like the ones 
presented in Zhang et al. (2009) and Beuther & Henning (2009) 
have recently started to throw light on scales where gravity is 
likely to dominate. 
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